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We explore the classical version of the mapping, due to Abanov and Wiegmann, of Calogero- 
Sutherland hydrodynamics onto the Benjamin-Ono equation "on the double." We illustrate the 
mapping by constructing the soliton solutions to the hydrodynamic equations, and show how 
certain subtleties arise from the need to include corrections to the naive replacement of singular 
sums by principal-part integrals. 
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I. INTRODUCTION 

The Calogero-Sutherland family of models consist of point particles moving on a 

line or circle and interacting with a repulsive inverse-square potential. Both the classical and 
quantum versions are completely integrable, and are the subject of an extensive literature 
3|. The models have application to one-dimensional electron systems 4j, and to the two- 
dimensional quantum Hall effect jlj]. 

It is possible to consider a hydrodynamic limit of the Calogero-Sutherland models in 
which the distribution and velocity of the particles are described by continuous fields p(x) 
and v(x) respectively. The equations of motion of these fields possess solitary wave solutions, 
and also periodic solutions that interpolate between small amplitude sound waves and large 
amplitude trains of solitary waves p,0|. The quantum version of the hydrodynamics provides 



an extension of the usual theory of bosonization of relativistic (i.e. linear dispersion) electron 
systems to systems where band-curvature effects become important js]. Because of the 
singular nature of the Calogero-Sutherland interaction, the hydrodynamic limit is rather 
more subtle than one might expect. Considerable insight into this limit has been provided 
by Abanov and Wiegmann who have shown 9| that it is equivalent to a "doubled" version 



of Benjamin-Ono dynamics 10|, llll]. The Benjamin-Ono equation, a member of an infinite 



hierarchy of integrable partial differential equations, was originally introduced to describe 
waves in stratified fluids. The connection has lead to a number of novel predictions for the 
evolution of one-dimensional electron gasses fl2 |. 

In this paper we will apply the tools of Q| to recover the classical soliton solutions found in 
6|, 13] and, in doing so, illustrate the structure of the Calogero-Sutherland Benjamin-Ono 
mapping. A secondary aim is to compare the origin of the hydrodynamic-limit subtleties in 
the classical model with their origin in the quantum system. In the quantized model they 
arise because of the need to convert operators that are hermitian with respect to one of the 
two natural inner products on the Calogero-Sutherland Hilbert space into operators that are 
hermitian with respect to the other 13j. In the classical model the subtleties arise because 
we need to include corrections to the naive limit of singular sums. 

In section[TTl we review the connection between the Calogo-Sutherland and Benjamin-Ono 
equations of motion. In section [TTT1 we consider an appealing, but overly naive, version of the 
hydrodynamic limit and reveal its failings. In section[lV]we illustrate how these shortcomings 
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are eliminated by including the first non-trivial order in an asymptotic expansion for the 
velocity field. In section |V] we consider the general mapping. An appendix contains proofs 
of some results used in the main text. 



II. CALOGERO-SUTHERLAND FROM THE BENJAMIN-ONO POLE ANSATZ 



The classical Benjamin-Ono equation [10, 111] is the nonlinear and and non-local partial 
differential equation 

u + u d x u = \X(d 2 xx u) H , (1) 
where /# denotes the Hilbert transform of /: 

p r°° i 

/hW = - -J(0d£. (2) 

If we introduce a Poisson bracket 

{u(x),u(x')} = 2\TrdJ(x-x'), (3) 

and Hamiltonian 

1 f°° { 1 A 1 

H *o = ^J_J- 6 u'- r (u x ) H ^dx, (4) 



then ([T]) can be written as ii(x,t) = {H BO , u ( x ,t)} 
infinitely many Poisson-commuting integrals of motion 



Following 



his Hamiltonian system possesses 



151 ]. we seek solutions of ([T]) as a sum of poles 



AT M 

The poles at a,j( t), j = 1, . . . , N, lie below the real axis while the poles at bj(t), j — 1, . . . , M, 
lie above it. In [15j, the numbers of aj and bj poles were set equal and bj(t) = a*(t). These 
conditions were imposed to ensure that u(x, t) was real. We will not make these assumptions, 
so our field u(x,t) is not necessarily real- valued. 
We insert the ansatz (15]) into (HI) and use 



x — a 



Ima < 



x — a ) H j +i 



(6) 



x — a 



, Ima > 
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to find 



N .,. M .,; "1 ( N . , M W iV . , M . , 



^(x-a fc ) 2 "^(x-6 fc ) 2 J ~ \ j-sx-aj~ j^x-bjj (x - a k f ~ (x - b k f 

-4ib^^ + y trA^\=°- ( 7 ) 

All terms with \/{x — a«) 3 and l/(x — fej) 3 cancel directly. The remaining terms can be 
simplified by exploiting the identity 

11 11 11 11,, 

+ 7Z -T7— v?> 8) 



(x - Ci) (x - c 2 ) 2 (x - c 2 ) (x - Ci) 2 (c 2 - Ci) (a; - c 2 ) 2 (ci - c 2 ) (s - Ci) 2 

to rearrange them as a sum of lf(x — a;) 2 's and l/(x — 6j) 2 's with x-independent coeffcients. 
Now the set of l/(x — cij) 2, s and l/(x — 6«) 2 's is linearly independent, and the vanishing of 
their individual coefficients requires 12] 



A " A 



A * A 

i6 ^ t £(^)-S(^)' (10) 

We have found iV + M evolution equations for iV + M variables, and so the pole ansatz is 
internally consistent. 

We now compute dj by differentiating and then using equations (jUJ), ({TO]) to eliminate 
the dj and fej's. After some labour involving repeated use of (jSJ) we find that 

_2A 2 

Remarkably, the bj(t) do not appear. The only role of the b/s in the aj-pole dynamics is 
that the complex parameters bj(0) determine the (complex) initial velocities dj(0) of the 
aj(t)-poles. Once these initial conditions are established, the poles evolve autonomously 
according to ffTTT) . 

The N equations (fTTI) are a complex version of the Calogero model equations of motion. 
As for the real-a^ case, they can be derived from the many-body Lagrangian 

N 

I . 

L 



Calogero ^ ^ ^ ^ , 



4 




We may site our initial poles so that if we place a pole at z, then we also place one at 
z + 27m for any integer n. The resulting 2tt periodicity will be preserved by the subsequent 
evolution. The sums 

= 5«*(f). («) 

00 I 1 

£ (7T2^ = 4 cosec2 (I) ■ (14) 

n=— 00 v ' 

then allow us to write the evolution equations in the form 

ia, = £ 5 «rt (^) - £ £ cot (^Z2i) , (15) 

^ 7 fe=l ^ 7 

fc;fc^j V 7 fc=l V 7 

and 

9 fl ^ A 



9 % 1 4 ^ sin (oj - a k )/2 



If we regard the a, as angles, then (fPTl) is the equation of motion arising from the Sutherland- 
model Lagrangian 

V - \2 



for particles on a circle. 



-^Sutherland , / ■ • , , / . o 

2 tr 8 ^ sm (°< - %)/ 2 



The authors of 



151 ] made 6, they wished u(x,t) to be real. Being interested 



primarily in Calogero- Sutherland models, we instead desire that the cii(t) be real. To arrange 
for this, we modify the definition of the Hilbert transform appearing in (jTJ. Following 
Abanov and Wiegmann 9] we define the T-contour Hilbert transform of u to be 

u r (z) d ^ -/ -L, f(C) df, (19) 

where T is a simple closed contour on which z lies. All that was needed to establish the 
a-pole autonomy was that the l/(x — a$)'s and l/(x — 6j)'s be eigenfunctions of the Hilbert 
transform with eigenvalues of opposite sign. Now if we take V to encircle the real axis in a 
clockwise sense (as shown in Figured]) then, for z on the contour, 



FIG. 1: The contour T surrounds the part of the real axis on which the ai poles are found. 



when the dj poles lie within V and the h\ poles lie outside. We can therefore let the eij(O) 
lie on the real axis and distribute the 6-poles in the remainder of the complex plane in such 
a manner that the initial dj's are real. Once the intitial eij's and dj's are real, the Calogero 
evolution ensures that the a$(t) remain on the real axis. Thus, a complex Benjamin-Ono 
field u(z,t) obeying 



on T can provide real-axis Calogero- Sutherland dynamics. 

III. SHEPHERD POLES AND CALOGERO DENSITY- WAVE SOLITONS 

In this section we will apply the pole ansatz to obtain solitary wave solutions to the 
hydrodynamic limit of the real-axis Calogero equation. In this limit we let the number N of 
Calogero particles become infinite, while keeping their density finite. We then replace the 
individual Oj(t) by a smooth particle-number density p(x,t), and the individual velocities 
dj(i) by a smooth velocity field v(x). We assume that p and v vary slowly on the scale of 
the inter-particle spacing. We are now naturally tempted to approximate the discrete a-pole 
sums in the evolution equations 



ii + u d z u = \\{&I z u)y 



(21) 




(22) 



by integrals to get 




A f°° A 

ih ti) = E i—r- / t— ( 23 ) 

Here, d(x, t) = t>(x, i) is the velocity of the pole at x. We begin by exploring the consequences 
of this approximation to the discrete sum. We will see that it is not quite consistent, and a 
small but significant correction is needed. 
Consider an initial density fluctuation 

p(x, 0) = po + Pi(x), where Pl (x) = (^) — 1 — . (24) 

Because 

'A 



7T J X 2 + A 2 



■d£ = l, (25) 



7T j e + A 2 

the Lorentzian distribution p\{x) corresponds to a local excess of one particle near x — 0. 
The contribution 

A ( A\ 1 Arc , T „ x ,^„ s 

^ = -^-T^2. (Ima; = 0) (26) 



L^t-x \tt J £ 2 + A 2 s x 2 + A 2 

of the density fluctuation to ia is real, and so tends to push a(x) off the real axis. Its effect 
can be countered, however, by placing a solitary 6-pole at at b = iA. We then have 

■ •/ «\ Aa; A 

ia(x,0) = - , ■ A2 + 



x 2 + A 2 x — iA 
\x X(x + iA) 



x 



+ A 2 x 2 + A 2 



iAX 



x 2 + A 2 ' 

and the a poles have a purely real initial velocity 



(27) 



^,0)=d(*,0) = ^^. (28) 



They therefore stay on the real axis. 

The motion of the b = iA pole is obtained from 

Apo 



A 

7T + ^vrpo 



6 + 

^ + ^7rp . (29) 
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Thus, the 5-pole velocity 

h = lA + Xnpo (30) 
is also purely real. Since we know that the a poles stay on the real axis, it must be that 

the a-pole density profile keeps abreast of the b pole as it moves parallel to the real axis. A 

constant-shape solitary wave of density 

p(x,t) = p + ( — ) 1 .x 2 . , 2 (31) 

therefore moves to the right at the speed tWiton = b. 

We could have put the b pole below the axis. In that case both the a-poles and pulse 
envelope will move to the left. In either case, the envelope velocity is always higher than 
the speed of sound t> SO und = A7rp , and is faster when the pulse envelope is tighter. 

Observe how the b pole acts as a shepherd: its real-part contribution serves to keep the 
a poles from wandering off the real axis. 

The a poles are distributed in such a manner that the b pole sees the effect of their 
enhanced density as a mirror image of itself lying below the axis. The contribution of 
this image to the b velocity parallels the interaction of a b and a = b* pair in the original 
real-valued Benjamin-Ono equation studied in [jjj] — the only difference being that in the 
present case we must include in b the constant velocity v SOU nd = A7rp induced by the uniform 
Po background. The condition that d(0) be real is a linear equation linking the fe-pole 
contribution to the density fluctuation 5p = (p — p ). Consequently the initial data for a 
chiral M-soliton can be established by placing poles above the axis at bj, j — 1, . . . , M. The 
associated Sp is then the sum of individual Lorentzians centered at Xj = Refej, and induces 
image poles below the axis at b*, j = 1, . . . , M. Thus, in our present approximation, the 
chiral Calogero hydrodynamics multisoliton coincides with the multisoliton solution of the 
conventional Benjamin-Ono equation. In particular, the real field 

^ t] = ^ { x^W) - ^W) } + X7Tp0 (32) 

obeys ((Tj), and possesses the physical interpretations 

u(x, t) = v + X-np = 2v + v sonnd = \-n(2p - p ). (33) 

This approximate mapping of the right-going Calogero density waves onto the Benjamin- 
Ono solitons is very appealing. Unfortunately there is a fly in the ointment: the equation of 
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continuity is not satisfied. For v and p both being functions of of x and t in the combination 
x — v so \i ton t, particle-number conservation requires that 



p + d x (pv) = d x {p{v - Ronton)} = 0. (34) 
Because p — > po an d v — > at large distance, the right-hand side of Q34D implies that 

pKoliton - U) = Po^soliton, (35) 

or, equivalently, 

v = Ronton- (36) 

P 

Now the a-pole velocity and density that we have found are linked by v(x) = \ir(p(x) — po), 
and since f so iiton is fixed while the p in the numerator varies with position, our present 
solution can obey fl36l) only approximately. 

IV. THE NON-LINEAR CORRECTION 

To get the continuity equation to hold exactly, we need to improve on the crude approx- 
imation 

£ > p (37) 

This naive continuum approximation would be legitimate (and fl36l) would hold exactly) were 
we to simultaneously let the background density po become infinite and take A — > in such 
a way that poA remains constant. We are not doing this, however. We wish to keep A fixed 
while po becomes large, but remains finite. For finite A, the poles immediately adjacent to 
Oj make a significant contribution to the sum, and their effect has to be carefully accounted 
for. 

An improved approximation to the pole sum is 

A f°° \ A 

^—p(U)dd--d x lnp(x)\ x=aj . (38) 

The d x In p correction is the first term in an asymptotic series that expands the difference 
between the sum and integral in local gradients of p. It arises because the particle at aj no 
longer lies midway between its neighbours at aj±i when the density varies with position - 
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but the symmetric cutoff in the principal-part integral tacitly assumes a midpoint location. 
We provide a derivation of this first correction term in the appendix. 

Once we know of the local correction to the principal-part integral, we realize that the 
shepherd poles need to make a corresponding additional real contribution to ia if they are 
still to prevent the aj poles from wandering off the real axis. The equations determining 
this extra contribution are non-linear, and so the simple pole superposition that enabled us 
to find the multisolition initial conditions in the previous section is not longer valid. We 
can still find some solutions, however. As an illustration, again consider the initial density 
profile 



P(x, 0) = Po + ( — 



7T J X 2 + A 2 



We now have have 



p f A A Xx X | A -2x 

P / 7 — rP(f,0)df- -a, hip = 



-oot-x 2 x 2 + A 2 2 I vr (x 2 + A 2 ) 2 _ , (A _ 

71 J X 2 + A 2 

Xx 



x 2 + A 2 + A/npo 
Xx 



x 2 + B 2 ' 



(39) 



where 



B=JA 2 + — . (40) 
ttPo 

We see that to keep the a's on the real axis, the shepherd 6-pole must be relocated to b = iB. 
The resultant a-pole motion is then governed by 

A _ A„ , A 



/°° A A 
■oo s x 



x — iB' 



Xx X 
+ 



x 2 + B 2 x-iB 1 
Xx ^X(x + iB) 



x 2 + B 2 x 2 + B 2 
iBX 



x 2 + B 2 ' 

The new, improved, v(x) is therefore 



(41) 



»w = (42) 
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From this, and (136|) . we get 



^soliton — v(x) 

P~ Po 



= PoAvr f^j , (43) 

which is independent of x as it should be — and was not previously. By squaring this last 
equation we also find that 

Anp = , (44) 

soliton sound 

which is the relation between soliton width and velocity obtained in |6j. 

Does this newly-computed soliton velocity coincide with that of the shepherd b pole? We 
set b = iB in 

ib = - X + iXnpo, (45) 
o + 1 A 

to find 

; A 



, • 'Vil-'- (46) 

Multiplying both sides by A(A + B) shows that b = iWiton, so the b pole does indeed track 
the density pulse — as it must if the pulse is to retain its shape. This consistency check 
also illustrates the fact that the d x In p correction does not affect the contribution of the 
a-pole sum away from the real axis. Once we are further away from the axis than the mean 
a-pole spacing, we no longer see the granularity of the a-pole "charge" distribution, and so 
the naive integral suffices for calculating the a-pole contribution to the 6-pole motion. 

We can also construct periodic solutions. If we take these to have period 27T, the resulting 
wave trains can be regarded as solitons of the trigonometric Sutherland model. Suppose, 
therefore, that the initial density profile is 



p(x, 0) = po + Yl 



n=— oo 



A\ 1 



7r / (x + 27m) 2 + A 2 

sinh A 
2tt J cosh A — cos x 



, 1 \ sinh A . , 



Then 



and 



-(A/2)sina; . , . 

-p(£,0)rf£= l u V , (48) 



£ — x ' cosh A — cos x ' 



. i T , A_ , — (A/2)sinx . . 

cosh A — cos x H 

2npo 
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This real contribution to ia can be cancelled by placing an infinite train of shepherd poles 
at z n = iB + 27m. From 

J ^ > ^ t( X ~ ^\ 1 /sinx + isinhS\ 

M— k» 1 (x-iB) + 27rn J ~ 2 C ° \ 2 / 2 \ cosh 5 - cos x J ' 

we see that I? should be chosen so that 

cosh I? = cosh A H sinh A. (51) 

By algebra that parallels the single-soliton case, we find that 

/ A \ sinh B 

v (x, = - — > 52 

v ' ; V 2 / cosh 5 - cosx V ; 

and that the soliton velocity is given by 

P i \ x f sinh £ \ . . 

^soliton = v(x) = poATl , . (53) 



and 



For a wave train of period A, these equations become 

v(x 0) = f^) ^H2nB/A) 

1 ' J V A J cosh(2 7 r J B/A) - cos(27rx/A)' 1 ' 

sinh(27r£?/A) 

Soliton = P0 A7F • r / I 55 ) 

sinh(27Ty4/A) 

They reduce to (1421) and (|43|) . repectively, when A becomes large. 

V. THE GENERAL CASE 

We now show that the local correction guarantees consistency with particle conservation 
for arbitrary initial data. 

To motivate the discussion, we begin by considering a simple model for the dynamics of a 
one-dimensional gas of spinless, unit-mass, fermions. The uniform non-interacting gas would 
have internal energy density e(p) = h 2 7T 2 p 3 /6. To include the effect of some interactions, we 
generalize this expression by introducing a parameter A so that e(p) = A 2 7r 2 p 3 /6. This A 
will later be identified with the parameter A appearing in the earlier sections. The simplest 
classical Galiliean-invariant Hamiltonian for the gas is then 

#fluid=y {i^ 2 +^-p 3 ] dx - ( 5e ) 
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Here v is the fluid velocity. We can write v = d x 9, where 9(x,t) is a phase field canonically 
conjugate to the density. The field 9 has Poisson bracket 

{9(x),p(x')} = 5(x-x), (57) 

leading to 

{v(x),p(x')} = d x S(x - x'), {v(x),v(x')} = {p(x),p(x')} = 0. (58) 
From p = {i?fl u id, p} we obtain the equation of continuity 

p + d x (pv)=0, (59) 
and from v = {H^^v}, we obtain Euler's equation 

(\ 2 2 2 \ 
■ (60) 

We can rearrange i?fl u id as 

and the equations of motion can similarly be massaged into Riemann form as 

d t (v + Xirp) + (v + \irp)d x (v + \irp) = 0, 

d t (v - Xttp) + (v- Xnp)d x (v - Avrp) = 0. (62) 

We therefore have non-communicating right-going and left-going Riemann invariants 
Ir,l{x) = (v±\irp) that are proportional to the chiral currents j'r ; l = \{p^ v I 'A7r) associated 
with the right and left Fermi points. The Riemann invariants have Poisson brackets 

{I R (x),I R (x')} = 2\irdJ(x-x'), 
{I L (x),I L (x')} = -2\ird x 6(x-x') : 

{I R (x),I h (x')} = 0. (63) 

Riemann's equations also show that this simple model contains the seeds of its own destruc- 
tion — the leading edge of any simple wave will inevitably steepen and break, making the 
fields unphysically multivalued. 

We now show how the Calogero gas mimics the Fermi gas, while regularizing the multi- 
valuedness. 
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Following [9J we decompose the u field in ^ as u(z,t) = u + (z,t) + U-(z,t), where 

M 

^t) = Y t J^ W y (65) 

The u±(z,t) are eigenf unctions of the T-contour Hilbert transform with eigenvalues ±z re- 
spectively. In the hydrodynamic limit, w_ becomes 

f°° A 

u-{z,t)=il p fat)j-—d£, (66) 



and has a discontinuity across the real axis: 

+ ie) — — ie) = 2\7rp(x). (67) 

Thus, 

w_(x ± ie) = \n(ipH ± p). (68) 
With the local correction included, the a-pole velocity v = a is 

f = %\npu + lnp + M + . (69) 

We can rearrange the last equation as 

u+(x, t)=v- i\np H - lnp. (70) 

For general p(x), v(x), we must take ( |70l) to be the definition of w + (x) on the real axis. We 
then define u+(z) to be the analytic continuation of u+(x) away from the axis. The analyt- 
ically continued u + will have no singularities within T, but will only be of the simple form 
(j65p for the restricted set of initial data that leads to pure soliton solutions. Nonetheless, 
as we show in the appendix, the only requirement for a-pole autonomy is that (w+)r = iu+. 
The total u = u + + u_ field therefore has limits above and below the real axis equal to 

u(x ± ie) — v ± A7rp — — In p. (71) 

These boundary values u(x ± ie) almost coincide with the Riemann invariants /r,l(^)- The 
only difference is a shift v — > v — i\d x lnp/2, and this shift does not affect the Poisson algebra 
in fl63|) . In the quantum theory, the invariance comes about because the shift is effected by 
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a conjugation 13], and conjugation does not alter c-number commutators. In the classical 



theory, we establish the invariance by observing that 

{6{x) - i\\\np{x),6{x') - p{x')} = 0, (72) 

and so, on differentiating by x and x', we find 

{v(x) -i\\d x hip(x),v(x') -i\\d x ,\np{x')} =0. (73) 

The shifted velocity therefore still Poisson-commutes with itself. Showing the invariance of 
the other brackets is straightforward. In particular, we find u(x + ie) has vanishing Poisson 
bracket with u(x — ie). Thus, the real-axis fluid-dynamics Poisson algebra (!58j) is equivalent 
to the natural extension of the Benjamin-Ono Poisson algebra ([3]) to the contour T. 
We can also reduce the T-contour version of the Benjamin-Ono Hamiltonian 

H r = -L I \-u 3 --u(d z u) r \dz (74) 



2A?r J r [6 4 

to a real-axis integral. Using the expressions flTTj) for u(x ± ie), and ur = i(u + — m_) we 
find, after some labour, that this integral is 

H hydro = J l-fnj 2 + —p 3 -—p(d x p) H + -^^Ux. (75) 

From the general theory, we know that (1741) leads to the Benjamin-Ono equation on T, and 
that that, in turn, leads to the Calogero equation on the real axis. Thus H hydTO , when used in 
conjunction with the real-axis fluid-dynamics Poisson algebra fl58l) . must be the Hamiltonian 
governing the real axis hydrodynamics. It therefore follows, from p = {if hydro , p} = —d x (pv), 
that the equation of continuity is exactly satisfied. 

-f^hydro is the continuum Hamiltonian of refs |f3, In those works it was obtained in- 
directly by taking h — > limit of the quantum continuum Hamiltonian that was already 
known from the collective field approach to matrix models [if]]. In the appendix we provide 
a more direct, and purely classical, derivation of the rather less-than-obvious expression for 
the internal energy part of -f^hydro : 

T ,r n lv- A 2 f°° ( A 2 7T 2 , A 2 7T ^ , A 2 (d x p) 2 1 , 

Although the boundary values u{x ± ie) of the u(z, t) field Poisson commute, the left- 
and right-going Riemann invariants are no longer dynamically decoupled. For example, an 
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inteplay between the two segments of the contour T is clearly seen when we use the w-field 
boundary values to compute the T-contour Hilbert transform of u at z = x — ie: 

( ■ \ P I u(z') J , 
Ut[x — ie) = — <p - — dz 



7r J r (x — ie) — z' 

1 f°° u(x' + ie) , , P f°° u(x' - ie) , , 

-ax / ax 



TT [x - ie) - (x' + ie) vr x - x' 



P f°° u{x' + ie) — u(x — ie) . 1 

-ax H — {nru(x + ie)\ 



IT J qq iC 3j IT 

P f°° 27cp{x r ) J , . , 

1 ' dx' + mix + ie) 



7T ./— oo % 



^npH + i (v + \7rp — i\\d x In p) 

Hu + -u-)\ z=x _ le . (77) 



(^) + m6{x - X>) (78) 



In the second line, the two integrals are the contributions from the upper and lower segments 
of T. The u(x + ie) appearing in the third line is the delta-function contribution from 

1 „ f 1 

x — x' — ie 

in the upper-segment integration. 

The interaction between the two branches of V also shows up when we seek chiral (i.e. uni- 
directional) waves. To have purely right-going motion, all the u + (z) singularities must lie 
in the upper half-plane. From the theory of Hilbert transforms, this condition is equivalent 
to demanding that 

(u+) H = iu+. (79) 
We can arrange for this eigenvalue condition to hold by imposing the chiral constraint 



0,Q,Q 



whence 



v = \ir(p-p Q ) - -X(d x lnp) H , (80) 

u + = Avr[(p - po) - ip H ] - -A [d x \np - i(d x hxp) H ] ■ (81) 

This last form of u + obeys (1791 because (f^H — —f f° r functions in L 2 (R). (Note that 
square-integrability requires us to subtract the background density from p in these equations. 
The Hilbert transform of a constant is zero, and (ph)h = —(p — Po)-) Under the chiral 
condition, the equation of continuity and the Euler equation reduce to a single equation 

p + d x [\Tc(p - po)p] = -d x {p{d x hxp) H \, (82) 
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This strongly non-linear equation reduces to the right-going Benjamin-Ono equation for the 
u = \n(2p — p ) held introduced in equations fl32|) and fl33l) if we linearize the dispersive 
term p(d x \iap) H — > (d x p)n, but in doing this we abandon strict number conservation. 



VI. DISCUSSION 



We have seen how the Benjamin-Ono equation on the "double" naturally generates the 
solitary-wave and wave-train solutions in the classical hydrodynamic limit of the Calogero- 
Sutherland model. We have also shown how this approach leads to an understanding of the 
the origin the subtle corrections to the naive continuous-fluid limit. We could have computed 
these corrections to higher order in the gradients of p. The small parameter in this expansion 
is, however, the local change in the inter-particle spacing divided by the spacing itself and 
should be small in the hydrodynamic limit. Keeping track of these higher-order terms would 
undo the advantages of the continuous-fluid approximation. We therefore retain only those 
that are required to maintain the internal consistency of the fluid mechanics model. 

It is interesting to compare the manner in which the classical machinery works with the 
quantum mechanical formalism of 0, [ljj. To make this comparison we should first note 
that the parameter A that appears in the quantum model is a dimensionless number. Our 
present classical parameter A has a scale, and is related to the quantum parameter by 

Aclassical ^Aq Uan t um . (^3) 

Thus, the free-fermion limit, where A quantum = 1, corresponds to A c i assica i = h, and any 
distinction between A quantU m = 1 and A quant um = is invisible in the classical H — > picture. 
What can be seen in the classical picture is the analaytic structure of the u(z, t) field and the 
resulting interplay of the left- and right-going currents Ir,l(^) — u(x ± ie) with the positive 
and negative parts u±(x) of the mode expansion. In particular, we see how these ingredients 
combine in a Poisson-bracket algebra that is the classical version of the quantum mechanical 
current algebra. 
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VIII. APPENDIX 

In this appendix we provide proofs of some assertions made in the main text. 



A. Singular sums 

Here we obtain the continuous-fluid approximations to the sums 

00 

51 = £ , (84) 

1 00 1 

52 = 2 £ (a m -a n y (85) 

n,m=— 00 v ' 

As usual, the singular a n = a and a n = a m terms are to be omitted in their respective 
sums. The essential tool is the Euler-Maclaurin expansion: 



00 poo 1 1 1 

E / w - y o fa) dv - - 12 /,(0) + 72d /,/,(0) + ' ' ' 



(86) 



This asymptotic expansion is valid for smooth functions f{y). The sums we are evaluating 
require singular /'s, however, and so a strategy of subtractions is needed. 

We regard the discrete numbers a n as the values at v = n e Z of a smooth monotonic- 
increasing function a(z/). We wish to form the sum 

Si = lim { V - - 1 - \ . (87) 

N ^ \z^ n a{n) _ a(0 ) J 



where the term with n = is to be omitted. We replace this sum by 

N 



r 1 

^ = lim V — -i 



a(0) a'(0)n J ' 



(88) 



where again the term with n = is to be omitted. Now S[ = Si because the subtracted 
sum is zero — all its terms cancel in pairs. Thus 



18 



where 

f{v) = a(u) - a(0) " ^ (90) 

has a smooth z/ — >■ limit: 

l^w - - < 91 > 

Keeping only the first correction in the Euler-Maclaurin series, we have 

J~ fly)dv ~ i/(0)+/(i) + /(2) + ..., 

f /(i/) &/ ~ J/(0) + /(-l) + /(-2) + • • • , (92) 

and so 

The integral in this last expression is convergent at v — since /(i/) is finite there. Because 
the Hilbert transform of a constant vanishes, we can, however, remove the counter-term 
that makes /(0) finite provided we cut off the v — > divergence with the principal-part 
prescription. On doing this, we obtain 

Si ~ P / — r^ dv + 



-00 
00 



= P / _^^__^ lnp(;r)L=ao (94) 



a(z/) - a(0) 2 [a'(0)] 5 
'-00 ? - a o 2 

In passing from the first line to the second we have made a smooth change of variables 
v 1 — ► £ — a(z/) (such changes of variables are legitimate in principal-part integrals) and 
defined the particle density p(x) in terms of a(v) by setting 

= £ = t = m < 95 > 

We have also used 

... . d 2 a 1 dp 1 dp 

tt » = = = (96) 

To evaluate the double sum S 2 , we begin by obtaining the continuum approximation to 

Ss= ^ (a(n)-a(O)) 2 (97) 

n=— oo v v ' \ / / 

where the n = term is to be omitted. Following our subtraction strategy, we use Euler's 
formula Y^=i l/™ 2 = vr 2 /6 to write 

1 7T 2 



V 



WfT=&{EW + ^» • (98) 
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where 



1 11 a "^ 1 

(Z/) ~ (a(u) - a(0)) 2 ~ KO)P^ + [^(Ojp^ ^ ' 



has been constructed to have a finite i/ — > limit: 



, r,/ \ 1 a"'(0) 3 a" 2 

m Fi/ = — r ,\ ; + - . 100 

i~o V ; 3 a' 3 4 a' 4 V 7 



We apply the Euler-Maclaurin formula as before to find that 

We can again omit the the 1/V term in F(y) at the expense of replacing the convergent 
integral by a principal-part integral. 
We now claim that 

pff ' - i 'UpflfWzfi ( 102 ) 

This result may be obtained by making a substitution v \— ► x = a(^) in the first term of 
the integrand on the left- hand- side, and a substitution v f— ► x = a(0) + a'(0)z/ in the second 
term. Now it is not immediately obvious that the use of two distinct changes of variable is 
legitimate. The integrals of the two terms do not exist separately - only the integral of their 
difference is convergent at v = 0. In order to be sure that no additional finite contribution in 
introduced by our manoeuvre, we must provide a common \u\ > e cutoff for the two separate 
integrals, and then keep track of the effect of the subsequent changes of variables on their 
integration limits. Because the two changes of variables agree to linear order near v = 0, we 
find that no such finite additions are induced. Our manoeuvre is indeed allowed. We also 
note that 

P r = t( p r f 1 *) = -*W*>»> = < 103 > 

J -oo Is % ) ^ V J —oo S & / 

Next we re-express F(0) as 

1 [a"(0)] 2 1 d ( a"{v) 



F(0) 



(104) 

v=0 



4 [a'(0)] 4 3rfz/ \[a'{u)f 
To complete the evaluation of S2 we should replace a by a n and sum over n. In the 

continuum aproximation we replace a(0) by a{v) and integrate over v. The total derivative 

in F will not contribute to this integral and can be discarded. After changing variables 

v — > x, we therefore find that 

S2 ~ Il{^ p3 ~l pidxp)H + l^} dX + '''' (1 ° 5) 
which immediately gives equation (176]) . 
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B. Pole-autonomy for a general « + . 



We here show that we can relax the condition that u + be a sum of simple poles, yet still 
have the aj obey the Calogero equation ffTTl) . 
Let 

N ., 

as before, but assume of u + (z, t) only that (u+)r = +iu + . Insert u = u_ + u + into 

u + ud z u=\d 2 zz u T . (107) 

The projections of the product u + U- onto the ±i eigenspaces of the Hilbert transform are 
respectively 

(«+«_)+ = - 

i=1 z a i 

iV 

«A 



(«+«-)- = Y] u + ( aj ). (10$ 



By taking the 2; derivative of (11081) . we we can project the cross terms u + d z u^ + u^d z u + 
appearing in (I107p into the ±z eigenspaces . From the coefficients of l/(z — cij) 2 in the — i 
eigenspace, we read off from (11071) that 

A 



+iu+(a j ). (109) 



, , , . Ofe — a j 

fc; fe^tj J 



In the +i eigenspace, we find that (11071) requires that 



N iX N i\ 
u + + u + d z u+ + V + V - {u+{a k ) - u+(z)) = i\\d 2 zz u + . (110) 



Now differentiate dj again to find that 

A 



Z („, _„.\2 ^3 - «fe) (^+( a i) + %^-U+| z = aj ) • (HI) 
k; k^j 



[a k - aj) 



If we momentarily forget all about u+, we know that we can assemble the remaining terms 
to obtain the Calogero equation 

2A 2 



Wi ( a i-°*) s ' 
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We therefore need to show that all terms in idj involving u+ drop out. These terms are 

7 — Z — ^2 ( n +( a i) -u + (a k ))+i [«+(%) + ^ — ~ — 9 zU+\ aj +u + d z u + \ a] J . (112) 
fc;Mj l a fc a j) V k-,k^j a i ak J 

Now consider the limit of ( 11 101) as z — > a,. There are potential singularities arsising from 
the /c = j terms in the sums, but on expanding 

u + (z) = u+(a,j) + {z- aj)d z u + \ a:j + \{z - a j ) 2 d 2 zz u+\ aj + 0[(z- c^) 3 ] , 

we find that the potentially-singular parts cancel among themselves, and, futhermore, the 
remaining finite parts of the j — k terms combine to cancel the iX^d zz u + term on the right 
hand side. The z = aj limit of f 1 11 j) thus equates to zero precisely the terms (11121) that we 
wish to disappear. 
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